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Abstract. A phenomenological model for the clustering of dark matter 
halos on the light-cone is presented. In particular, an empirical prescrip- 
tion for the scale-, mass-, and time-dependence of halo biasing is described 
in detail. A comparison of the model predictions against the light-cone 
output from the Hubble Volume A^-body simulation indicates that the 
present model is fairly accurate for scale above ~ 5/i~^Mpc. Then I ar- 
gue that the practical limitation in applying this model comes from the 
fact that we have not yet fully understood what are clusters of galaxies, 
especially at high redshifts. This point of view may turn out to be too 
pessimistic after all, but should be kept in mind in attempting precision 
cosmology with clusters of galaxies. 



1. Introduction 

The standard picture of cosmological structure formation suggests that any vis- 
ible object forms in a gravitational potential of dark matter halos. Therefore, 
a detailed description of dark halo clustering is the most basic step toward un- 
derstanding the clustering of visible objects in the universe. For this purpose, 
many theoretical models for halo clustering have been developed and then tested 
against extensive numerical simulations. 

First, I will describe our most recent theoretical model for clustering of dark 
matter halos (Hamana, Yoshida, Suto & Evrard 2001b). In particular, we focus 
on their high-redshift clustering where the past light-cone effect is important. 
Then I will show that our model predictions are in good agreement with the 
result from a light-cone output of the Hubble volume simulation (Evrard et 
al. 2001). Finally I will discuss a fundamental difficulty in relating the halo 
model to clusters of galaxies. My conclusion is that we already have a reliable 
empirical model for the halo clustering, but that we need to understand what are 
the clusters of galaxies, especially at high redshifts, before attempting precision 
cosmology with clusters of galaxies. 

2. Theoretical model for two-point correlation functions of dark mat- 
ter halo on the light-cone 

As emphasized by Suto et al.(1999), for instance, observations of high-redshift 
objects are carried out only through the past light-cone defined at z = 0, and 
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the corresponding theoretical modeUng should properly take account of a variety 
of physical effects which are briefly summarized below. 



2.1. Nonlinear gravitational evolution of dark matter density fluctu- 
ations 

Assuming the cold dark matter (CDM) paradigm, the linear power spectrum of 
the mass density fluctuations is computed by solving the Boltzmann equation 
for systems of CDM, baryons, photons and (usually massless for simplicity) 
neutrinos. The resulting spectrum -Pii^icar(^' ^^^^ space is specified by a 
set of cosmological parameters including the density parameter Qq, the baryon 
density parameter fif,, and the Hubble constant h in units of 100/km/s/Mpc, 
and the cosmological constant Aq. Then one can obtain its nonlinear counterpart 
in real space, P^{k, z) by adopting a fitting formula of Peacock Sz Doods (1996). 



2.2. Empirical model of halo biasing 

The most important ingredient in describing the clustering of halos is their bi- 
asing properties. The mass-dependent halo bias model was developed by Mo 
& White (1996) on the basis of the extended Press-Schechter theory. Subse- 
quently Jing (1998) and Sheth &; Tormen (1999) improved their bias model so 
as to more accurately reproduce the mass-dependence of bias computed from 
A^-body simulation results. We construct an improved halo bias model of the 
two-point statistics which reproduces the scale-dependence of the Taruya & Suto 
(2000) bias correcting the mass-dcpcndcnt but scale-independent bias of Sheth 
&; Tormen (1999) on linear scales as follows: 

6haio(M,i?,z) = 6sT(M,z)[1.0 + 6sT(M,^)(7M(i?,2)]°•'^ (1) 

for R > 2Rvij-{M, z), where RviriM,z) is the virial radius of the halo of mass 
M at z. In order to incorporate the halo exclusion cfi"cct approximately, we set 
bhaioiM, R,z) = for < 2Ryij-(M, z). In the above expressions, (Jm{R-,z) is 
the mass variance smoothed over the top-hat radius R = (^M / AirpaY^^ , Po is 
the mean density, 5c{z) = 3{12TTy/^ /20D{z), D{z) is the linear growth rate of 
mass fluctuations, and v = [5c{z)/(Jm{R,z = 0)]^. 



2.3. Redshift-space distortion 

In linear theory of gravitational evolution of fluctuations, any density fluctua- 
tions induce the corresponding peculiar velocity field, which results in the sys- 
tematic distortion of the pattern of distribution of objects in redshift space 
(Kaiser 1987). In addition, virialized nonlinear objects have an isotropic and 
large velocity dispersion. This finger- of- God effect significantly suppresses the 
observed amplitude of correlation on small scales. With those effects, the non- 
linear power spectrum in redshift space is given as 



P^{k,II,z) = P„^(fc,z)[l +/3halo/^^]^£'vel[fcMf^halo], 



(3) 
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where £)vei is the Fourier transform of the pairwise pecuhar velocity distribution 
function (e.g., Magira et al. 2000), /x is the direction cosine in fe-space, (Xhaio is 
the one-dimensional pair-wise velocity dispersion of halos, and 

_ 1 d\nD{z) 
Phalo = -7 Tj ■ [V 

bhaio d\nz 

While both cihaio and 6haio depend on the halo mass M, separation R, and z in 
reality, we neglect their scale-dependence in computing the redshift distortion, 
and adopt the halo number-weighted averages: 

/ 2D^ {z)a^ (M, z = 0)nj (M, z)dM 

c7Lio{>M,z) = ^ , (5) 

/ nj{M,z)dM 
Jm 

/•oo 

/ bsT{M,z)nj{M,z)dM 

bhaio(>M,z) = , (6) 

/ nj{M,z)dM 
Jm 

where wc adopt the halo mass function nj{M, z) fitted by Jenkins et al. (2001). 
The value of (T^{M, z = 0), the halo center-of-mass velocity dispersion at z = 0, 
is modeled following Yoshida, Sheth & Diaferio (2001). Then our empirical halo 
bias model can be applied to the two-point correlation function of halos at z in 
redshift space as 

ehaio(M, R, z) = 6Lio(M, R, z) P^k, ^)^^- (7) 
2.4. Cosmological light-cone effect 

All cosmological observations are carried out on a light-cone, the null hypersur- 
face of an observer at z = 0, and not on any constant-time hypersurface. Thus 
clustering amplitude and shape of objects should naturally evolve even within 
the survey volume of a given observational catalogue. Unless restricting the ob- 
jects at a narrow bin of z at the expense of the statistical significance, the proper 
understanding of the data requires a theoretical model to take account of the 
average over the light cone (Matsubara, Suto, & Szapudi 1997; Mataresse et al. 
1997; Moscardini et al. 1998; Nakamura, Matsubara, & Suto 1998; Yamamoto 
Sz Suto 1999; Suto et al. 1999). According to the present prescription, the two- 
point correlation function of halos on the light-cone is computed by averaging 
over the appropriate halo number density and the comoving volume element 
between the survey range Zmin < z < Zr, 



'-max 

j dM j dz-^ n]{M,z)Ui^{M,R,z) 

^b2i0{>M,R) = , (8) 

/ dM dz^n]{M,z) 
Jm Jz^^„ dz 



where dVc/dz is the comoving volume element per unit solid angle. While the 
above expression assumes a mass-limited sample for simplicity, any observational 
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selection function can be included in the present model fairly straightforwardly 
(Hamana, Colombi & Suto 2001a) once the the relation between the luminosity 
of the visible objects and the mass of the hosting dark matter halos is specified. 




Figure 1. Two-point correlation functions of halos on the light-cone; 
simulation results (symbols; open circles and filled triangles for real 
and redshift spaces, respectively) and our predictions (solid and dotted 
lines for real and redshift spaces, respectively). The error bars denote 
the standard deviation computed from 200 random re-samplings of the 
bootstrap method. The amplitudes of ^^'^ for Mhaio > 5.0x lO^^/i'^M© 
are increased by an order of magnitude for clarity. 

Figure 1 compares our model predictions with the clustering of simulated 
halos from "light-cone output" of the Hubble Volume ACDM simulation (Evrard 
et al. 2001) with Qb = 0.04, J^cdm = 0.26, as = 0.9, Qa = 0.7 and h = 0.7. For 
the dark matter correlation functions, our model reproduces the simulation data 
almost perfectly at i? > 3/i^^Mpc (see also Hamana et al. 2001a). This scale 
corresponds to the mean particle separation of this particular simulation, and 
thus the current simulation systematically underestimates the real clustering 
below this scale especially at z > 0.5. Our model and simulation data also 
show quite good agreement for dark halos at scales larger than 5/i~^Mpc. Below 
that scale, they start to deviate slightly in a complicated fashion depending 
on the mass of halo and the redshift range. This discrepancy may be ascribed 
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to both the numerical limitations of the current simulations and our rather 
simplified model for the halo biasing. Nevertheless the clustering of clusters on 
scales below 5h~^Mpc is difficult to determine observationally anyway, and our 
model predictions differ from the simulation data only by ~ 20 percent at most. 
Therefore we conclude that in practice our empirical model provides a successful 
description of halo clustering on the light-cone. 



4. Dark matter halos vs. galaxy clusters 

With the successful empirical model of halo clustering, the next natural question 
is how to apply it in describing real galaxy clusters. In fact, in my opinion the 
main obstacle for that purpose is the lack of the universal definition of clusters. 
Let me give some examples that I can easily think of (see Fig.2). 



Abell clusters 

selected by eyes 

from the photographic plates 

galaxies satisfying 



m < 2 



Press-Schechter halos 

nonlinear spherical collapse 
A„. = 18 7t' 



SZ clusters 

inverse Compton of CMB 
AIsz~neTeRci 



Halos in N-body simulations 

friend-of-friend method 
linking length 

= 0.2 (mean particle separation) 




X-ray clusters 

thermal bremsstrahlung 

1/2t 



n"T""R 



Figure 2. Dark halos - galaxy clusters connection. There are a vari- 
ety of practical definitions of dark matter halos and clusters of galax- 
ies. They are deinitely related, but the one-to-one correspondence is 
unlikely and nothing but a working hypothesis. 



(i) Press-Schechter halos: almost all theoretical studies adopt the definition 

of dark matter halos according to the nonlinear spherical model. This is 
characterized by the mean overdensity of IStt^ (in the case of the Einstein 
- de Sitter universe. The corresponding expressions in other cosmological 
models can be also derived.). Combining this definition with the Press- 
Schechter theory, the mass function of the dark halos can be computed 
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analytically. This makes it fairly straightforward to compare the predic- 
tions in this model with observations, and therefore this definition has 
been widely studied in cosmology. 

(ii) Halos identified from N-body simulations: in reality the gravitation- 

ally bound objects in the universe quite often show significant departure 
from the spherical symmetry. Such non-sphcrical effects can be directly 
explored with N-body simulations. Even in this methodology, the identifi- 
cation of dark halos from the simulated particle distribution is somewhat 
arbitrary. A most conventional method is the friend-of-friend algorithm. 
In this algorithm, the linking length is the only adjustable parameter to 
controle the resulting halo sample. Its value is usually set to be 0.2 times 
the mean particle separation in the whole simulation which qualitatively 
corresponds to the overdensity of ISvr^ as described above. 

(iii) Abell clusters: until recently most cosmological studies on galaxy clus- 
ters have been based on the Abell catalogue. While this is a really amazing 
set of cluster samples, the eye-selection criteria applied on the Palomar 
plates are far from objective and cannot be compared with the above def- 
initions in a quantitative sense. 

(iv) X-ray clusters: the X-ray selection of clusters significantly improves the 
reliability of the resulting catalogue due to the increased signal-to-noise 
and moreover removes the projection contamination compared with the 
optical selection. Nevertheless the quantitative comparison with halos de- 
fined according to (i) or (ii) requires the knowledge of gas density profile 
especially in the central part which fairly dominates the total X-ray emis- 
sion. 

(v) SZ clusters: the SZ cluster survey is especially important in probing the 

high-z universe. In this case, however, the signal is more sensitive to the 
temperature profile in clusters than the X-ray selection, and thus one needs 
additional information/models for temperature in order to compare with 
the X-ray/simulation results. 

The above consideration raises the importance to examine the systematic com- 
parison among the resulting cluster/halo samples selected differently. In reality, 
this is a difficult and time-consuming task, and one might argue that we do not 
have to worry about such details at this point. Such an optimistic point of view 
may turn out to be reasonably right after all. Nevertheless it is still impor- 
tant, at present, to keep in mind that this simplistic assumption of "dark halos 
= galaxy clusters" may produce a systematic effect in the detailed comparison 
between observational data and theoretical models. 

5. Conclusions 

I have presented a phenomenological model for clustering of dark matter halos 
on the light-cone by combining and improving several existing theoretical mod- 
els (Hamana et al. 2001b). One of the most straightforward and important 
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applications of the current model is to predict and compare the clustering of 
X-ray selected clusters. In doing so, however, the one-to-one correspondence 
between dark halos and observed clusters should be critically examined at some 
point. This assumption is a reasonable working hypothesis, but we need more 
quantitative justification or modification to move on to precision cosmology with 
clusters. 

I am afraid that this problem has not been considered seriously simply 
because the agreement between model predictions and available observations 
seems already satisfactory. In fact, since current viable cosmological models 
are specified by a set of many adjustable par anient crs, the agreement does not 
necessarily justify the underlying assumption. Thus it is dangerous to stop 
doubting the unjustified assumption because of the (apparent) success. I hope 
to examine these issues in future. 
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